Comparison of implicit and explicit AUSM‐family schemes for compressible multiphase flows

This paper makes the first attempt of extending implicit AUSM‐family schemes to multiphase flow simulations. Water faucet, air–water shock tube and oscillating manometer problems are used as benchmark tests with the generic four‐equation two‐fluid model. For solving the equations implicitly, Newton's method along with a sparse matrix solver (UMFPACK solver) is employed, and the numerical Jacobian matrix is calculated. Comparison between implicit and explicit AUSM‐family schemes is presented, indicating that similarly accurate results are obtained with both schemes. Furthermore, the water faucet problem is solved using both staggered and collocated grids. This investigation helps integrate high‐resolution schemes into staggered‐grid‐based computational algorithms. The influence of the interface pressure correction on the simulation results is also examined. Results show that the interfacial pressure correction introduces numerical dissipation. However, this dissipation cannot eliminate the overshoots because of the incompatibility of numerical discretization of the conservative and non‐conservative terms in the governing equations. The comparison of CPU time between implicit and explicit schemes is also studied, indicating that the implicit scheme is capable of improving the computational efficiency over its explicit counterpart. Copyright © 2014 John Wiley & Sons, Ltd.     

POD enriched boundary models and their optimal stabilisation

Modes obtained using the Proper Orthogonal Decomposition are used as boundary enrichment functions within a variational multiscale method for the stochastically forced Burgers equation. Initially, large increases in accuracy are obtained using the enrichment functions without stabilisation terms. Then, optimal coefficients for the stabilisation parameter τ of the unresolved scale model are calibrated using a goal‐oriented model‐constrained optimisation technique, resulting in further improvements. As both the determination of the enrichment functions and the optimisation of the coefficients requires high‐accuracy reference data, a scaling procedure is introduced to allow their use over range of conditions. Numerical experiments confirm that the scaling procedure is effective. Copyright © 2014 John Wiley & Sons, Ltd.     

Molecular structure and properties of click hydrogels with controlled dangling end defect

In this study, controlled amount of dangling ends is introduced to the two series of poly(ethylene glycol)‐based hydrogel networks with three and four crosslinking functionality by using click chemistry. The structure of the gels with regulated defect percentage is confirmed by comparing the results of low‐field NMR characterization and Monte Carlo simulation. The mechanical properties of these gels were characterized by tensile stress–strain behaviors of the gels, and the results are analyzed by Gent model and Mooney–Rivlin model. The shear modulus of the swollen gels is found to be dependent on the functionality of the network, and decreases with the defect percentage. Furthermore, the value of shear modulus well obeys the Phantom model for all the gels with varied percentage of the defects. The maximum extension ratio, obtained from the fitting of Gent model, is also found to be dependent on the functionality of the network, and does not change with the defect percentage, except at very high defect percentage. The value of the maximum extension ratio is between that predicted from Phantom model and the Affine model. This indicates that at the large deformation, the fluctuation of the crosslinking points is suppressed for some extend but still exists. Polymer volume fractions at various defect percentages obtained from prediction of Flory–Rehner model are found to be in well agreement with the swelling experiment. All these results indicate that click chemistry is a powerful method to regulate the network structure and mechanical properties of the gels. © 2016 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 1227–1236     

Sorption of methanol,dimethyl carbonate,methyl acetate,and acetone vapors in CTA and PTMSP: General findings from the GAB Analysis

Sorption of vapors of four organic compounds in two glassy polymers, cellulose triacetate (CTA) and poly[(trimethylsilyl)propyne] (PTMSP), has been reported and analyzed in terms of Guggenheim‐Anderson‐De Boer (GAB) model. These two structurally and physicochemically different glassy polymers both independently showed that one sorption site was formed by about three monomeric units. This finding held true for vapors of all characterized compounds; that is, for methanol, for its derivatives dimethyl carbonate and methyl acetate, and for acetone. The “rule of three” might thus also be applicable to other sorbates and glassy polymers. Further, an original modification of the GAB model for the sorption of alcohols in PTMSP was derived and successfully tested. Overall, the analyses of the sorption isotherms, heats of sorption and diffusion coefficients supported the view that the sorption of vapors in glassy polymers has adsorptive nature. © 2015 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2016 , 54, 561–569     

超音速等离子体炬的磁流体动力学数值研究

针对设计的喉径2mm、工作电流为100A的拉瓦尔喷嘴,在二维轴对称模型的基础上,对超音速等离子体炬中的流动及其外部射流进行了数值模拟。通过在阳极喷嘴内部采用基于磁矢量势的磁流体动力学模型,避免了对磁感应强度的复杂积分计算,得到了喷嘴内部多场耦合的结果及外部射流的流动状态,分析了喷嘴内部电磁场对等离子体的加速作用及射流发展过程。结果显示,等离子体经历了亚音速→跨音速→超音速的发展过程,最终获得2.3 Ma的超音速射流。研究结果为超音速等离子体炬的工业应用提供了理论基础。     

A review of computational modeling and deep brain stimulation: applications to Parkinson's disease

Biophysical computational models are complementary to experiments and theories, providing powerful tools for the study of neurological diseases. The focus of this review is the dynamic modeling and control strategies of Parkinson's disease (PD). In previous studies, the development of parkinsonian network dynamics modeling has made great progress. Modeling mainly focuses on the cortex-thalamus-basal ganglia (CTBG) circuit and its sub-circuits, which helps to explore the dynamic behavior of the parkinsonian network, such as synchronization. Deep brain stimulation (DBS) is an effective strategy for the treatment of PD. At present, many studies are based on the side effects of the DBS. However, the translation from modeling results to clinical disease mitigation therapy still faces huge challenges. Here, we introduce the progress of DBS improvement. Its specific purpose is to develop novel DBS treatment methods, optimize the treatment effect of DBS for each patient, and focus on the study in closed-loop DBS. Our goal is to review the inspiration and insights gained by combining the system theory with these computational models to analyze neurodynamics and optimize DBS treatment.     

干式变压器有限元仿真模型的电磁和振动分析

随着电力需求的逐年增长,干式变压器的数量也在不断增加。干式变压器在运行时存在着振动和噪声的问题,为了对干式变压器振动的规律与特点进行研究,本文建立了干式变压器本体振动的有限元仿真模型,通过电磁分析获得相应的磁场分布,然后利用结构动力学分析得到其本体振动的相关规律。通过对处于运行状态的变压器振动数据进行实测分析,得到变压器振动的特征频率,然后对仿真结果进行对比分析,可以发现振动幅度与频率之间存在的关系。本文的研究结果可对干式变压器的减振降噪研究提供参考。     

Acoustoelastic effect simulation by time–space finite element formulation based on quadratic interpolation of the acceleration

Acoustoelastic effect describes the change of ultrasound velocity due to the initial stress. Its simulation involves a numerical analysis of nonlinear elastodynamics and requires high accuracy in the time domain. A time–space finite element formulation, derived from the quadratic interpolation of the acceleration within a time segment, is proposed for an accurate simulation of the acoustoelastic effect in the present study. Ten different integration schemes are generated based on this formulation and nine of them are found to be conditionally stable. Among the nine stable schemes, one is found to obtain a spectral radius of one when the normalized step ratio is less than 5.477, indicating no numerical dissipation or numerical divergence. Compared with integration schemes from previous studies, this integration scheme demonstrates better performance in calculation accuracy and energy conservation. A two-stage approach, namely the static stage and the dynamic stage, has been employed in the simulation of the acoustoelastic effect. The former stage is adopted to obtain the initial stress and the latter stage, where the proposed integration scheme is implemented, is adopted to simulate the ultrasound propagation in an initial stress state. The simulation results of the dynamic stage show that the ultrasound velocity increases in a compression stress state and decreases in a tension stress state for aluminum alloy, which is in good agreement with previous experimental studies. Together with the simulation result of the static stage, it is conjectured that the acoustoelastic effect results from the stress-dependent elastic modulus.     

双模随机交错晶场中最近邻弱交换相互作用对重入现象的影响

利用有效场理论研究了双模随机交错晶场中混合自旋 Blume-Capel 模型纳米管系统的重入现象, 发现了系统的重入现象与晶场取值概率、 晶场强度和外壳层与内壳层格点间最近邻交换相互作用的关系. 结果表明: 取值概率、 交换相互作用、 晶场强度和温度等诸多因素相互竞争, 使系统表现出丰富的磁化现象: 正( 负) 晶场较弱时, 系统只发生二级相变; 随着正( 负) 晶场增强, 系统的二级相变消失, 呈现一级相变; 一定条件下, 系统会出现重入现象.